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Abstract 

We demonstrate, by an explicit one-loop calculation, that at leading twist the 
nonperturbative effects in i? — i? — 77 and B radiative decays 

are contained in a common multiplicative factor {Ab{E^), where E^^ is the energy 
of the photon). We argue that this result holds also at higher orders. Ratios of 
the amplitudes for these processes do not depend on scales below the mass of the 
5- meson (Mb), and can be calculated as perturbative series in asiMs). 



1 Introduction 



In this letter we make the observation that, at leading power in the mass of the i?-meson 
(Ms), the nonperturbative QCD effects in the radiative 5-decays, B — > -yiiye, B — > 'y£~^i~ 
and S — > 77 are universal. By this we mean that, for the same photon energy E-y, the 
amplitudes for the processes B — > ^ii^g and B — >■ ^£'^£~ are proportional to each other. 
For = Mb/2 they are also proportional to the amplitude for B '-fj. The constants of 
proportionality contain kinematic factors, CKM-matrix elements and a series in ^^(Mb) 
which is therefore calculable in perturbation theory. Thus, in spite of the different weak- 
decay mechanisms for the three processes, the strong interaction effects are common at 
scales below Mb- 

In rcf. [1] we showed that the two form factors (Fj^y) for the B ^ii'e decay are given, 
up to one-loop order, by the factorization formula: 

Fa,v{E^) = J dk+^l{k+- fip) TCk+, E^, fip) , (1) 



with 

■^SCET( \ fsQuMB \ 1 



2^/2 E^ J k. 



1 H Kt{k+, Ej] jiF) 

47r 



(2) 



where the momentum of the photon is in the — direction, C^^^'^ {j-ii. ) is a Wilson coeffi- 
cient relating the weak b ^ u current to operators of the Soft-CoUinear Effective Theory 
(SCET) [2], is the S-meson's leptonic decay constant and — —2/3 is the charge of 



the iZ-quark ^. hf is a factorization scale and is conveniently chosen to be 0{^y Mb Aqcd)- 
The S-meson's light-cone distribution amplitude, is the component of the non-local 
matrix element [4, 5] 

<^ap{h) = J dz.e^'^^- {0\M^)MbamH)l^^^^^, , (3) 

with the hadronic initial state H — B, which contributes to the amphtudes at leading twist. 
In eq. (3), a, (5 are spinor labels and [z,Q] denotes a path-ordered exponential. Finally, 



,r ^ ^ . o2V2E^k+ TT^ 
Kt{k+, E^- Pip) = log' - - 1 ■ 

Hp D 



The above formulae hold for photon energies satisfying E^ ^ Aqcd- 

The nonperturbative contribution to the form factors is contained in ^^{k^) which 
is convoluted with T as in eq. (1). The observation we make here is that precisely the 
same convolution over k+ appears also in the amplitudes for B — > ji'^i" and S — > 77 
decays. The A;+-independent prefactor in curly brackets in eq. (2) is process-dependent but 

^Although we have chosen to use the SCET formaUsm to sum the large logarithms, this could equiva- 
lently be achieved using ear Her approaches, in particular the "Wilson-Line" formalism of ref. [3]. 
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Figure 1: Lowest-order diagram at leading twist contributing to the process bu — > 7X, 
where X = iug, or 7 and is represented by the dashed hne. The thick hne repre- 

sents the 6-quark and the thin hne a hght quark. The grey circle represents the operator 
responsible for the h ^ u transition. 

perturbative. Moreover, although the coefficient C^^^"^ {^p) in eq. (2) is specific to the 
decay B — > 7^//^, the behaviour with (for < itib) is common to all coefficients C'^'cbt 
appearing in the three radiative decays. Thus ratios of amplitudes for different processes 
but with the same only depend on the Cf^^^ at fip = Mb, and are hence calculable 
as perturbative series in q;s(Mb). The formulae for the decay rates are presented explicitly 
in sec. 3. 

Additional motivation for this study is provided by the need to understand nonpertur- 
bative QCD effects in two-body exclusive nonleptonic decays, for which there is a wealth 
of experimental data, particularly from the 5-factories. The demonstration of the factor- 
ization of long-distance effects at leading twist for these processes [6, 7, 8] provides the 
theoretical framework for detailed phenomenological analyses. However, our ignorance of 

(and the lack of a similar framework for higher-twist contributions) limits the precision 
with which information about the CKM matrix elements and CP-violation can be deter- 
mined from the experimental measurements of branching ratios and asymmetries. It is 
therefore important to understand whether there are relations, similar to those presented 
here, also between different nonleptonic decay amplitudes (or between contributions to the 
amplitudes) which reduce or eliminate the need for a detailed knowledge of Finally we 
repeat that measuring the photon's energy distribution in B ^ 'yii^i decays may in future 
be the best way of obtaining information about [1]. 

In the next section we sketch the evaluation of the amplitudes for i? — > 77 and B 
'y£~^i~ decays at one-loop order and the resummation of the large logarithms. We focus on 
the differences with the corresponding calculation of the B —>■ jiui decay amplitude, which 
was described in detail in ref. [1]. The reader who is only interested in the implications 
of our results can turn immediately to section 3 where we present the expressions for the 
decay rates and discuss their significance. 



2 Amplitudes up to One-Loop Order 

In ref. [1] we have described the evaluation of the one-loop contribution to the hard- 
scattering kernel for the decay B 'yii'e in detail. The corresponding calculation for the 
other radiative decays is very similar, so here we only sketch the main points and present 
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Figure 2: One-loop leading-twist diagrams for the process bu — X7. 
the results. 

We wish to write the form factors for each of the three radiative processes in the generic 
factorized form given in eq. (1), where the hard-scattering kernel T depends on the process, 
but, as its name suggests, does not depend on scales below fip and can be calculated in 
perturbation theory. In order to determine T we are free to choose any appropriate and 
convenient external state, and we take a w-antiquark with momentum k and a 6-quark with 
momentum p-k. The leading- twist contributions are those in which the photon is emitted 
from the light quark, and we will see that the diagrams which need to be evaluated are 
those in figs. 1 and 2. 

The evaluation of the one-loop graphs in fig. 2(a),(d),(e) and (f) is common to all 
the processes and is described in ref. [1]. In ref. [1] we also show that the leading- twist 
contribution from the box diagram in fig. 2(c) comes from the soft region of phase-space 
and is absorbed into the one-loop component of the wave-function, giving no contribution 
to T. This is also true for S — » 77 and B — ^ '~f£'^i~ decays. Thus in the following we 
present the results for the diagrams in fig. 1 and 2(b) for 5 — >■ 77 and B decays 
and combine them with those in ref. [1] for the other graphs (the graphs are evaluated in 
the Feynman gauge). 



2.1 The decay B — ^ 77 

The effective Hamiltonian for 5° — > 77 decays is [9, 10] 

n = tXl'^^^m,Cf{f,n) [qYYbR]{d,A, - d.A,) , (5) 

where Br = 1/2(1 -|- 75)6, q = d or s and A^"^^ = Vt^th- Cf^ is the Wilson coefficient at the 
scale /iij (it will generally be convenient to take /x/j = M^). The superscript eff denotes the 
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fact that Cj^ contains terms proportional to the Wilson coefficients Ci-q (corresponding 
to AB = 1 four-quark operators) and Cg (corresponding to the chromomagnetic operator). 
The contributions included in Cj^ are short- distance loop corrections to Oi_6 and Og 
which result in an effective vertex proportional to the operator O7 at tree-level (i.e. to the 
operator in eq. (5)). 

We now briefly illustrate how the leading- twist contributions proportional to Ci_6 and 
Cg are absorbed into C^^. It is not surprising that the diagram in fig 3(a) has the structure 
of fig. 1, where the grey circle represents the electromagnetic operator O7, and its contribu- 
tion can therefore be included into C^^. Indeed this is conventionally done [9]. In fig. 3(b), 
(c) and (d) we draw, for illustration, three diagrams at 0{as) contributing to C^^. The key 
point to note is that the loop integrations are dominated by the short-distance region of 
phase space, and therefore the contributions to C^^ are calculable in perturbation theory. 
In fact Cj^ at 0{as) can be deduced from the calculations in refs. [11, 12] (for the penguin 
operators O^^q, for which the Wilson coefficients are small, the contribution to Cj^ is only 
known at 0{a^)). In the basis and notation of refs [11, 12] 



7 ^ Q ^ Q zLtt * ^ ' ^ ' 

1,2,8 



3 9 3 9 47r 



where and Fg are given in eq. (56) of ref. [11] and eq. (82) of ref. [12] respectively. 
There are other contributions proportional to Ci_6 and Cg which cannot be absorbed into 
C^^ . An example is drawn in fig. 3(e). However, power-counting arguments show that 
these contributions are suppressed by at least one power of Mb- The contributions to 
the hard-scattering kernel from the 0(0;^) diagrams drawn in fig. 2, where the grey circle 
represents the tree-level insertion of O7, are evaluated explicitly in this letter. 
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The amplitude (A) for the decay S — > 77 can be written in terms of invariant form 
factors A±: 

A{B 77) = ^^{2^- el el)+A^ [2^{q^ ■ el){q, ■ el) - zM^e^ • 6^) ] }. (7) 

Following ref. [1], we choose the external state H = bq to determine the contribution 
to the hard-scattering kernel. We denote the contribution from the tree diagram in fig. 1, 
where the photon with momentum qi is emitted from the effective vertex, to the {'y^\'H\bq) 
matrix element by F^^^ (0 denotes tree-level or "zero loops"). Similarly ^2°^ is the contri- 
bution with the photon with momentum q2 emitted from the effective vertex, and is equal 
to Fi^^ with the obvious replacements 1^2 and + ^ — in the expressions below. 

Evaluating the diagram we find for the matrix element at tree level 

^''^ ^ ^''^^7^"^'^'^ ■ VfS;^'''^^^ m2[ivm + 75)] U^p - k), (8) 

where u and v are the wave functions of the b and q quarks with spin labels S and s. 
Writing this amplitude as a convolution 

(0) _ f rf/c+ feg- (0)^(0) /q^ 
(with ^''^ defined in eq. (3) with H = bq), the contribution to the hard-scattering kernel is 

Adding the contribution from T2^''^, we obtain for the lowest-order contribution to the two 
form factors: 

A+ = A_^ jdk+^^{k+;i,F)T'-''\k+;i,R,i,F)^xi'^Cf{i,R)^, (11) 
r(o) - C^^(a ^^^'^ ^ - ^ I dk im 



Note that is also the nonperturbative quantity which is needed in evaluating the con- 
tribution from tree-level hard-spectator interactions to exclusive nonleptonic i?-decay am- 
plitudes [6, 7, 8]. 

At one-loop order, the contribution of the diagram in fig. 2(b) to the {'y^\'H\bu) matrix 
element is 



^(l)wk ^ ^(0) O^sCf 



_V^-21og^-vr^ 
^ Mb Mr 



(13) 
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In order to determine the one- loop contribution to T (T"^"^^), we need to subtract the 
corresponding contribution to $''^(^) (g) T^'^\ the convolution of the one- loop component of 
with r(°): 



(0) 



4n 



4 2 2A;2 
+ - log- 



9 I '■^a 2 



-2 log' 



(14) 



where we use dimensional regularization, working in 4 — cu dimensions. We recall that the 
subscript 1 denotes the contribution from the graphs with the photon of momentum qi 
emitted from the effective vertex. 

The contribution to the hard-scattering kernel from the diagram in fig. 2(b) with qi ^ q2 
is the same, since the resulting convolution over is equal by symmetry to the one over 

from fig. 2(b). We therefore present our answers below in terms of a convolution with 

as the (dummy) integration variable, and our expressions for the hard-scattering kernel 
correspond to the sum of the two contributions. 

The results in eqs. (13) and (14), together with those in ref. [1] for the remaining 
diagrams in fig. 2 can be combined to give the one-loop expression for the hard-scattering 
kernel: 



1 + 



An 



log' 



y^k^ 
M 



+ 4 log 



^/2k^ 
M 



1 ^ 
log — 



M 



M 



+2 log^ h71og — 



2 log — 

jlF 



(16) 



where the explicit expression for T^^^ is given in eq. (12). 

We now use the SCET formalism to resum the large logarithms following the steps in 
sec. 5 of ref. [1]. The contribution of the tensor operator to the {'^^\T-L\hu) matrix element 
(i.e. the diagrams in figs. 1 and 2(b) plus half of the diagrams in figs. 2(d) and 2(e)) is 



An 



M, 



2 log log 



V2k+ 7 , ^/2k^ 

log 

Mb 2 ^ Mb 



— - — TT 



OLsCl 

4n 



01 1 i 

2 log log 

IJ'R 



V2k. 
~M^ 



7, 
log 

2 ^ 



V2k^ 



M, 



n 



B 



(16) 



In the SCET formalism, only one operator, denoted by Og in ref. [2], contributes to the 
form factors. From eq. (66) of ref. [1] and the neighbouring discussion, we can rewrite F"^^ 
in terms of the SCET Wilson coefficient C|^'^^^(E^; hr, hf) as 



:^Jw 



CI^''^{Mb/2;^r,Mb) X 



(0) 



47r 



log 



2 V2k^ 



7, ^/2k^ 
- log 

2 ^ Mb 



7 

+ 2- 



11^^ 
"T2" 



+Ff ) I 1 + 



An 



log 




(17) 
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Comparing eqs. (16) and (17) we find 



2 47r 



21og 6 - — 

piR 12 



(18) 



which agrees with eq. (33) of ref. [2] for hji = Mb- 

Finally, we can use the renormalization group equation for Cg'-^^^, derived in ref. [2], to 
obtain its value at any factorization scale /ip. We conclude that the form factors — A_ 
are given by the generic form in eq. (1) with 



(19) 

Thus the hard-scattering kernel has exactly the same dependence on A;+ as that for B jii/^ 
decays in eq. (2). 

2.2 The decay B ji+£- 

The procedure for the evaluation of the B decay amplitude follows exactly the 

same steps. We will assume here that = 2E^/Mb satisfies (1 — x^) 3> Aqcd/-^s- This 
allows us to neglect the diagrams in which the real photon is emitted from the weak vertex 
and the virtual one is radiated from the hght quark. The effective Hamiltonian for this 
process is then [9, 13, 14]: 

7rv2 V Qi 

(20) 

where bi.B. = (1 Tl5)/2b. Again Cf^ contains short-distance, leading-twist contributions 
from four-quark and chromomagnetic operators: 

Cf^Cg + Y{ql)-^Y.C^F^''^ (21) 

1,2,8 

where Y{qf) has contributions from Oi-e and its explicit expression is given in cq. (10) 
of ref. [12] and Fi^2 ^-nd Fg^^ are given in eqs. (54) and (55) of ref. [11] and eq. (83) of 
ref. [12] respectively. We follow the authors of ref. [14] in neglecting the long-distance 
effects associated with cc resonances. For the purposes of this paper, we start with eq. (20) 
and evaluate the matrix elements of the operators on the right-hand side. The four form 
factors for this process are defined by: 

2(7(?, e*)\q^,bL\B{p)) = ee^.,,e* Wi^y(^7) + ' <l) ' ' ^1]Fa{E,) (22) 

and 

2{^{q,e*)\qia^,{p-qybR\B{p)) = -eMBe^,,,e* ^^^^(^7) 

-ieMB[el{v ■ q) - q,{v ■ e*)]F^{E^) . (23) 
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Since the corrections to the vector and axial currents have aheady been discussed in ref. [1], 

we consider here only the form factors Fj- and F^. Evaluating the diagram in fig. 1 we find 
that at lowest order the contribution to the matrix element 2{t^ l~ '^\qi(T nvqihR\hq) is 



ir(o) 



2(g2)-A;H 

and the hard-scattering amplitude is 

7^(0) ^ 



r^2[7M,^l](l + 75)] 



Pa 



The form factors at lowest order are thus given by 

Fr = F^^ ^^^^^^ 



(24) 



(25) 



(26) 



where = 1/3 is the electric charge of q. 

At one- loop order the contribution from the diagram of fig. 2(b) to the matrix element 

2{e+e-^\qia^,q-(bR\bq) is 



-log^--21og--21og-log--log^--2Li2 1 - - -tt 

y y y y y \ x/ 



, (27) 



where x = M^, y = 2MbEj and z = 2q2-k. The corresponding contribution to (g)T^°^ 
is given by the r.h.s. of eq. (14) with Ff '' replaced by F^^^ given in (24). 

Combining the above results with the remaining diagrams in fig. 2 from ref. [1] we find 
that the hard-scattering kernel is given by 



T 



l + log' - + 2 log - log 4 - 2 log ^ - 2 log^ - 

An \ X X Hp y y 

1,92; ^-r . y\ , H% 3 , X TT^ „\ 



(28) 



We now perform the resummation of the large logarithms as for the B '-ff decays 
above. By comparing our results for the contribution of the tensor operator to the matrix 
element F'^^' with the results for the same transition in the SCET, we obtain the expression 
of Cg'^'^'^{Ej; Hp — Mb) (it is given explicitly in eq. (33) below). For /jr — Mb our 
result agrees with eq. (33) of ref. [2] . 

In this way we find that the form factors Ft and F^ are given by the generic form in 
eq. (1) with 



T{k+,E^;HR,Hp) 



IbMb 1 ^SCET 

6V2E^ ~k 



[E^] Hr, Hf) 



■7 ^+ 



1 + ^^^Kt{k+,E^;HF) 



(29) 



This completes the demonstration that at one-loop and leading-twist order the hard- 
scattering kernel has the same form for all three radiative decays. 
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3 Discussion and Conclusions 

As a consequence of the universality of the nonperturbative effects in radiative i?-decays 
the rates can be related simply using perturbation theory. The decay rates are given by 
the following expressions: 



r{B, ^ 77) 

dr{Bq 7e+e-) 
dK 



144vr3 



1 



A|(Mb/2) ' 



(30) 



17287r4 



Al{E. 



X 



7V 



(jeS (jSCET _|_ 



9 /^cfT 

1 

1 — X 



7 



+ Cio C3 



SCET I 2 



547r2 



.r-, (1 — x^) 



(31) 
(32) 



where = 2E^/Mb- C^^^'^ in eqs. (30) -(32) denote the matching coefficients relating 
the QCD and SCET operators evaluated a,t fip — Mb and are calculable in perturbation 
theory, 



riSCET 

^9 



and 



(jSCET 



1 + 



= 1 + 



as{MB)CF 
An 



a,{MB)CF 
An 



7J-2 

log - 2 log^ x^ + 2 log x^ - 2Li2 {1 - Xj) - 6 - — 



(33) 



21og^a;^ - 2Li2 (1 - x^) + 



n 



logx^-e . (34) 
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jiR is the scale where the Wilson coefficient functions of the operators in the effective 
Hamiltonian, Cf^. Cf^ and Cio, are evaluated, and for convenience wc set fiR = Mb- The 
only nonperturbative parameter in eqs. (30) -(32) is Ab{E^), which is the generalization 
of Xb: 



1 



1 



= — P exp[—S{EY, ij,f)] / dk. 



1 + 



An 



, 22V2E^k+ n^ , 
log ^2^- y-l 



The resummation of large (Sudakov) logarithms is obtained through the evolution of SCET 
operators from Mb down to /if (this evolution is common to all the operators), leading to 
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the exponential factor in Ab{E^) 

Mr) 



-S{E^;I1f) 
Mr) 

Mr) 



as{mi,) 



PI 



+ Mr) 



1 + log r 

r 



1 — r + r log r — 



(36) 
(37) 



log^ r 



/9o 



— 2 log 



logr — 



2 

2CfB 



[r - 1 - logr] , 



(38) 



where r = a^i^ip) / ois{mb). Po = 11Ca/3 - 2iV//3 and /?i = 34Cy3 - lQCANf/3 - 2CFNf. 

In this letter we do not perform a detailed phenomenological analysis. For illustration 
however, in fig. 4 we plot the ratio 



Rq{E^) - 



dT{Bq 'ye+e-)/dE^ 



(39) 



dr{B+ ^e+u)/dE^ ' 

for q = d,s and for the following values for the Wilson coefficients: Cf^ — —0.390, Cg^ = 
4.182 and Cio = —4.234 [9, 13, 14] (a more detailed analysis would include the effects 
of l^(gi) and the Oia.^ corrections to these coefficients, but these are sufficiently small 



to be neglected in this presentation). We take |A 



7.8 X 10"^ |A 



4.1 X 10-^ 



iKbl =3.6x 10-^ 



Ra X 



190 MeV and /b, 

2 



230 MeV. In other words we plot 
' and 




(40) 



(41) 



The two curves are proportional to each other, differing by the values of the leptonic decay 
constants and CKM matrix elements. They do not depend however, on the model for the 
distribution amplitude. As approaches 1 the virtual photon approaches its mass shell 
and this is the reason for the rise of the curves in this region. We recall however that our 
expressions are valid in the region 1 — Aqcd/^s (and x^ ^ Aqcd/Ms). 

Knowledge of A^, which we do not have at present, is required to determine each of the 
branching ratios separately. Taking Ab{Mb/2) = 350 MeV, we estimate that the B —>■ jj 
branching ratios are approximately 

77) ~ 10"^ , 



B{Bd ^ 77) ~ 3 X 10" 



and 



B(B.. 



(42) 



of the same order of magnitude as in ref. [10]. For the remaining decays we estimate: 

B{B+ ^£+u) = O(10~^) , (43) 
B{Bs 7^+r; x^ < 0.8) = 0(10"^) , (44) 



B(B, 



<0.8) = 0(10"") -0(10"^°), 



(45) 
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Figure 4: Ratios of decay rates Rd (solid line xlO~^) and Rg (dashed line xlO~^) as 
functions of = 2E^/Mb- 



where we have imposed the cut-off < 0.8 in order to satisfy the constraint {1 — x^) ^ 
^qcd/Mb- Assuming that there is no very large enhancement as ^ 1, the result for 
B{Bs '~fi^i~) is in broad agreement with the estimates in ref. [13]. This result is an 
order of magnitude smaller than that of ref. [14] . The discrepancy is partly due to the fact 
that we do not reproduce the formula for the ratio B{Bs 'ji~^i~)/B{Bs —>■ 77) used in 
ref. [14], partly to a different choice of input parameters, partly due to the fact that the 
authors of [14] integrate up to the kinematic limit for x^ and partly to the fact that we 
include higher-order QCD corrections. 

The explicit calculations described in the previous section, were performed at one-loop 
order. Nevertheless, as we now explain, we anticipate that the general features will survive 
also at higher orders. For all the decays the lowest-order diagram of fig. 1 gives a factor of 
l/k+ which leads to the factor of 1/^+ in the hard-scattering kernels. At one-loop order all 
the diagrams in fig. 2 except for 2(c) have a collinear light-quark propagator external to the 
loop giving a factor of In order for the diagram in fig. 2(c) to contribute at leading 

twist we must recover the factor of 1/A;+ from the loop integration. As explained in ref. [1], 
such a contribution comes from the soft region of phase space and is absorbed into the 
wave function - it cancels when ^^^^ (g) T^^^ is subtracted, and therefore does not contribute 
to the hard-scattering kernel. Thus the one-loop contribution to the hard-scattering kernel 
has two components: 

• the correction to the electromagnetic vertex - figs. 2(a),(e) and half of (d)), 

• the correction to the weak vertex - figs. 2(b), (f) and half of (d)). 

The power-counting arguments leading to this structure appear to be sufficiently general 



11 



for us to conjecture that it is true to all orders. The corrections to the electromagnetic 
vertex arc then clearly common to all three decay processes. Even though the corrections 
to the weak vertex are not common, the differences only arise at scales of 0{Mb). For soft 
and coUinear gluons, the Dirac structure of the heavy and light quarks is indeed such that 
the result is independent of the form of the weak operator, and of the 7-matrix (matrices) 
in particular. Although there are two large scales in each of these decay processes, Mb and 
^/ MbAqcb, the differences arc only due to physics at the scale M b and can be expressed 
in terms of a calculable perturbation series in ^^(M^). Recent studies of factorization at 
higher orders of perturbation theory include a demonstration of the validity of eq.(l) for 
B ■jii'i decays [15] and a study of the transitions from soft to coUinear quarks in the 
SCET [16]. 

The processes studied in this letter have been rare radiative decays. For other pro- 
cesses, such as the two-body exclusive decays of B-mesons into two light mesons [6, 8], the 
hard-scattering kernels will in general be different. Nevertheless, as a result of the inde- 
pendence of soft and coUinear QCD effects from the structure of the weak operator, it is 
to be expected that there are analogous relations between contributions to different decay 
amplitudes, and we arc currently investigating the structure and scope of these relations. 

We end by repeating that in this letter we have only considered the leading-twist 
contributions to radiative 5-decays. For a detailed and precise phenomenology of exclusive 
S-decays the extension of the factorization formaUsm to the 0(Aqcd/Mb) corrections will 
be necessary (see, for example, refs. [17, 18] for recent contributions towards this). 
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